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We derive recursively the perturbation series for the ground-state energy of the D-dimensional 
anharmonic oscillator and resum it using variational perturbation theory (VPT). From the exponen- 
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orders. The calculation effort is much smaller than in the standard field-theoretic approach based 
, _ on the Hubbard-Stratonovich transformation. 
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i The properties of nontrivial physical systems can only be calculated via efficient approximation schemes. Most 
easily accessible are perturbation expansions, but they are usually divergent and need rcsummation. To this end, 
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I. INTRODUCTION 



ft' 

a variational approach was developed by Feynman and Kleinert which has been systematically extended to an 
efficient nonperturbative approximation scheme called variational perturbation theory (VPT) It allows the 

conversion of divergent weak-coupling into convergent strong-coupling expansions and has been applied successfully 
in various fields, such as quantum mechanics, quantum statistics, condensed matter physics, and critical phenomena. 
In fact, the most accurate critical exponents come from this theory as has been verified by recent satellite 
\ experiments 0. 

m 

The convergence of VPT has been analyzed up to very high orders for the ground-state energy of the one-dimensional 
"f^ , anharmonic oscillator 

Oh" 1 

V(x) = -u; 2 x 2 +gx A (1) 

c : 

and was found to be exponentially fast 0, 0. This surprisingresult has been confirmed later by studying other 
physical systems and was proven to hold on general grounds Furthermore, the exponential convergence seems 

to be uniform with respect to other system parameters. In this manner, the variational resummation of perturbation 
series yields ap proximations which are generically reasonable for all temperature s [ HI . ITU H^ . space and time 
coordinates 0,0,0,0, magnetic field strengths 0], coupling constants [l8L Il9l l20| . etc. 



In this paper, we show that the exponential convergence of VPT is uniform with respect to the space dimension D. 
To this end, we consider the D-dimensional generalization of the anharmonic oscillator l(T|l. i.e. 

V(x) = ic 2 x 2 + 5 (* 2 ) 2 (2) 

with x = (x\, . . . , xd), and determine its ground-state energy as a function of the coupling constant g. In Section II, 
we derive the corresponding weak-coupling series by evaluating connected vacuum diagrams. In Section III, we show 
how this perturbation series can be obtained more efficiently by means of Bender- Wu-like recursion relations |2l| . 
In Section IV, we resum the weak-coupling series by applying VPT and examine the resulting convergence, which is 
again exponentially fast and improves uniformly with increasing dimension D. In Section V, we show that the latter 
observation is not surprising, since the ground-state energy of the anharmonic oscillator (J2J can be determined with 
the help of a systematic large-D expansion. In Section VI, we apply VPT to extract the large-D expansion to higher 
orders, which have so far been inaccessible to other methods. 
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II. PERTURBATION THEORY 

The perturbation series for the ground-state energy of the anharmonic oscillator can be calculated from connected 
vacuum diagrams. Up to the second order in the coupling constant g, the ground-state energy is given by the Feynman 
diagrams 



S = --lim4l-00 

2 p^oo (3 \ 8 ^- A - ; 



1 



1G 



ooo 



(3) 



with the propagator 



!2i e -oj\T a -T b \ 
2u 



(4) 



and the vertices 



SiiSki + SikSji + SuS-jk }■ I dr a . 



(5) 



The connected vacuum diagrams © are derived most easily by an efficient graphical recursion method j2^| . Evaluating 
these Feynman diagrams produces the following analytic expression for the ground-state energy: 



E 



Dlo D{D + 2)g D(2D 2 + 9D + I0)g 2 



4j 2 



8lo 5 



(6) 



Only low perturbation orders are accessible by this procedure. If we want to study higher orders, we better use 
Bender- Wu-like recursion relations [2l| . 

III. BENDER- WU-LIKE RECURSION RELATIONS 

The potential J5J of the Z?-dimensional anharmonic oscillator is rotationally symmetric. Hence, the ground-state 
wave function depends only on the distance x — |x|. We solve the corresponding Schrodinger eigenvalue equation 



1 / d 2 D — 1 d\ 1 2 2 4 

2{d^ + —rd^ + 2 UJX +9X 



ip(x) = Etp(x) 



(7) 



as follows. We write the wave function ip(x) as 



^0*0 = {-) ex P 



lx 2 + 4>(x) 



(8) 



with the abbreviation x = Xy/uj, and expand the exponent in powers of the dimensionless coupling constant g = g/u 
by using 



= £<M*)s* 



(9) 



fc=i 



The 4>k{x) are expanded in powers of the rescaled coordinate x: 



fc+i 



(10) 



For the ground-state energy we make the ansatz 



E = LU [ T + E efc 5 fc 



(11) 



fc=i 
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k 




1 


D(D + 2)/4 


2 


-D(2£> 2 +9£> + 10)/8 


3 


D(8D 3 + 59D 2 + 146D + 120)/16 


4 


-D(168D 4 + 1773D 3 + 7144D 2 + 12960D + 8840)/128 


5 


L»(1024D 5 + 14325D 4 + 82222D 3 + 241464D 2 + 360736L> + 216960)/256 



TABLE I: Expansion coefficients for the ground-state energy of the anharmonic oscillator up to the fifth order. 



Inserting JSJ-^U into (0, we obtain to first order 



D + 2 



1 



4 ■' "2 --J> £ i 



£>(£> + 2) 



For fc > 2, we find the following recursion relation for the expansion coefficients in H10|) 

fc— 1 m 



(to +!)(£> + 2m) (k) 
2m Cm4 



n(m + 1 - n) (l) (fc _rj 



1=1 n=l 



(k) 

with = for to > k + 1. The expansion coefficients of the ground-state energy follow from 

n (fc) 
e fe = -Dc\ . 



(12) 



(13) 



(14) 



Table [I] shows the resulting coefficients e& up to the fifth order. For D = 1, they reduce to the well-known one- 
dimensional results I2II. 



IV. VARIATIONAL RESUMMATION 

Now we consider the strong-coupling limit of the perturbation series JHJ. Rescaling the coordinate according to 
x — > xg -1 / 6 , the Schrodinger equation becomes 



1 / a 2 D-ia 

2 \<9x 2 x dx 



1 



,- 2 / 3 „, 2 -r- 2 



+ x* 



rl>(x) = g-^ 3 Eip{x) 



Expanding the wave function and the energy in powers of the coupling constant yields 

tf>(x) = M%) + M%) g~ 2/3 + M%) g~ 4/3 + ■ ■ ■ , 

E = Lug 1 ^(bo + b 1 g- 2 ^ + b 2 r 4/5 + ---) ■ 



(15) 

(16) 
(17) 



By considering (|15fl in the limit g — > 00, we find that the leading strong-coupling coefficient bo equals the ground-state 
energy associated with the Hamilton operator 



Precise numerical values for this ground-state energy for different dimensions D are listed in Table ITT1 
The weak-coupling series (|ll|l is of the form 



(18) 



' n N 1 

T + !>(£) 



(19) 



k=l 



4 



b (D = 2) 


1.4771497535779972(31) 


b (D = 3) 


2.3936440164822970(37) 


b (D = 10) 


10.758265165443755(69) 



TABLE II: Numerical results for the leading strong-coupling coefficient bo for the ground-state energy of Q18I . 



b (D = 2) 


1.477149753577994356(33) 


b (D = 3) 


2.3936440164823030895(77) 


b (D = 10) 


10.758265165443797408091(18) 



TABLE III: High-precision VPT results from the 80th order for the leading strong-coupling coefficient bo in the ground-state 
energy 0. 



The alternating signs and fast growing coefficients in Table [I] suggest that (|19|l is a divergent Borel series which is 
resummable by VPT QE||1,|^. To this end, an artificial frequency parameter ft is introduced in the perturbation 
series most easily by Kleinert's square-root trick 



u CLy/l + gr, (20) 

with 

r = ^~- < 21 > 

One replaces the frequency uj in the weak-coupling series l|19|) according to (|20fl and re-expands the resulting expression 
in powers of g up to the order g . Afterwards, the parameter r is replaced according to l|21|) . This procedure has the 
effect that the power series of order N for the ground-state energy becomes dependent on the variational parameter 
ft: 



k=a i=o ^ . \- 



(22) 



The influence of ft is then optimized according to the principle of minimal sensitivity |23|: the ground-state energy 
to A^th order is approximated by 

£W=£M( S)W) !lW), (23) 
where ft^) denotes that value of the variational parameter for which E^ N ' (g, u>, ft) has an extremum or a turning point. 

As an example, consider the weak-coupling series (I19fl to first order: 

m Duj D(D + 2) , . 

Inserting l|2U|l . re-expanding in g to first order, and taking into account (|21|1 . we find 

4 ' 4ft + 4ft 2 

Extremizing this and going to large coupling constants, we obtain the strong-coupling behavior of the variational 
parameter 



m , , Dn Duj 2 D(D + 2) , , 

EW (g, u, ft) = — + -— + ^ ] g . (25) 



ft« = oj 5 1/3 (ft^ + ftl 1} .g- 2/3 + 4 1} g- 4/3 + ■■■), (26) 

with 

ft« = {/2(D + 2) , ftW = , Uj2 - , ftW = f s , .... (27) 

V V j ' 1 3^/2(^ + 2) ' 2 108(£> + 2)' V ' 
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Inserting the result l|26[) . (|27|l into l|25[) . we obtain the strong-coupling series (|17f) with the first-order coefficients 

#J = 2^0 + 2), M - DUj2 bP = .... (28) 

8 V 1 1 4^2(L> + 2) 2 48(1? + 2)' 1 ; 

For D = 1, this reduces to earlier VPT results in Refs. 0,0. 

In order to determine the strong-coupling coefficient bo in l|17|) to higher orders, we proceed as follows. We observe 
that the strong-coupling behavior of the variational parameter is for any order N of the form 

nw = ^ 1/3 (4 N) + n^g- 2/a + 4 N) g~ 4/3 + ■■■) (29) 

which corresponds to l|26|) . Inserting l|29|) into l|22|) . the leading strong-coupling coefficient turns out to be given by 

= ee ( (i " r )/2 ) (-i)^^) 1 -" - (30) 

where the inner sum can be further simplified by using |24| 



m / x 

/— n \ / 



-l) m (' a (31) 



Thus, the leading strong-coupling coefficient reduces to 

N 



111 



k=0 ^ ' 

Determining the optimized Og W ' for which b N \fl N ^) has an extremum or a turning point then leads to the 
approximation b (Qq ) of the leading strong-coupling coefficient 6o- 

By carrying the expansion to high orders, VPT yields approximations whose relative deviation from the exact value 
vanishes exponentially 0.[J. In our case, we have for large N 



\b ( N \D)-bo(D)\ ^ 
TTFn ~ ex P 



A{D) - B(D)N 1/3 , (33) 



where the exponent 1/3 is determined by the structure of the strong-coupling series l|17|) . Due to the exponential 
convergence of VPT, it turns out that the accuracy of our numerical results for the leading strong-coupling coefficients 
bo in Table[n]is not sufficient for a useful examination of the convergence behavior of VPT in high orders. Therefore, 
we use our results from the 80th VPT order as a more precise approximation for bo. Table ITTT1 summarizes our 
high-precision VPT results, where in each case the uncertainty of bo has been estimated by examining the deviation 
from the result of the previous order. For D = 2 and D = 10, the VPT results lie within the error margins of the 
numerical results. However, this is not the case for D = 3, where the VPT result lies just outside of the corresponding 
numerical error margins. We attribute this discrepancy to an overly optimistic error estimate for the numerical result. 
The precision of the results shown in Table II I II improves with increasing dimension, which already indicates that the 
calculation converges faster in higher dimensions. Figure ^ shows the convergence of the VPT results for the three 
different cases. Fitting the data to straight lines yields for the parameters A(D) and B(D) in l|33|) : 

A(D = 2) = 5.98(72) , B{D = 2) = 9.89(23) , (34) 
A(D = 3) = 7.43(48), B(D = 3) = 10.67(15) , (35) 
A(D = 10) = 11.89(63), B(D = 10) = 13.33(20) . (36) 

Thus, we find that the convergence of the VPT result indeed improves uniformly with increasing dimension D. 
This tendency does not come as a surprise, since for D — > oo the ground-state energy of an oscillator with quartic 
anharmonicity can be determined exactly as we will sec in the subsequent section. 



1 8 1 / 3 27 1 ' 3 64 1 / 3 

FIG. 1: Logarithm of the relative deviation of the VPT result for the leading strong-coupling coefficient bo(D) for D = 2 
(circles), D = 3 (diamonds), and D = 10 (triangles) plotted as a function of the cubic root of the perturbation order up to 
70th order. The dashed lines represent least squares fits of the data to straight lines. 



V. LARGE-D EXPANSION 



Now, we elaborate the systematic large- D expansion for the ground-state energy of the anharmonic oscillator @ 
based on standard field-theoretic methods (see, for instance, Refs. [3l 12a l26| ) . 



A. Effective Potential 



We start with the path integral representation of the quantum-statistical partition function at finite temperature T 

Z= <j)Vxe- A M , (37) 

where the Euclidean action reads 

rf> fi - 1 - - - .1 

(38) 



A[x] = j'dr |ix 2 (r) + ^V(r) + g [x 2 (r)f j 



The paths are periodic in the imaginary time t with period (3 = 1/T. Applying a Hubbard-Stratonovich transforma- 
tion, 



Da exp • 



dr 



-a 2 (T) + 2ix 2 (T)a(T) 
9 



exp < -g 



/ dr [x 2 (7 
Jo 



the x(r)-path integral l|37|l . (|38|l becomes harmonic and leads to 



Z = (bVae 



-DAM 



where we have introduced the Euclidean action 



- 4 -M = 7T / dr a 2 (r) + ~Tr In 
Dg Jo 2 



d 2 

— —r+u; 2 + Ai<j{T) 



(39) 



(40) 



(41) 



The remaining path integral (|40|l . (|41|1 over ct(t) is then performed in the limit of large D by regardin g th e modified 
coupling constant g = Dg as being independent of D. and by applying the background method pi. l27l l28| . Thus, we 
take into account order by order the effect of the fluctuations 8<j(t) = <t(t) — cto of the paths <j(t) from the background 
(To- This determines the effective potential 



in the form of the loop expansion 



y eff (ao) = £v$(a ) : 
/=( 

where the term of loop order I turns out to be of order D 1 



(42) 



(43) 



1=0 



B. Loop Orders I — and / = 1 



The leading term is the tree-level and follows from evaluating l|41|l for the background ctq 



r(0), 



with the auxiliary frequency 

The one-loop correction is given by 

where the operator 
turns out to be 



= D 


M + 


— In ( sinh 






. 9 


P V 




Q = 


v^ 2 


+ 4zcr • 





Kff (^o) = ^TrlnG- 1 (r 1 ,r 2 ) 



2/? 



(5cr(Ti)5cr(T 2 ) 



ct(t)=<t 



G" 1 (n , r 2 ) = - 6 (n - r 2 ) + 8 G 2 , (n , r 2 ) . 
9 

Here, the correlation function Gq(ti,t 2 ) has the spectral representation 

°^ e -iwm(n-T 2 ) 
Gu(n,r 2 )= ^ - , 

m— — oo v 7 

with the Matsubara frequencies cj m = 27rm//3. Inserting 1)49)1 in 1148)1 yields 



1 oo 

G- 1 (r 1 ,r 2 ) = - £ G ™ 



with the coefficients 



Gm — ~ + ^ 



-1 e — iw m (n— t 2 ) 



m— — oo 



m' — — cxd 



(w 2 ,+ft 2 )(w 2 



For the quantum-mechanical ground state energy to be calculated we only need the zero-temperature limit /3 
of the above quantum statistical expressions. Thus, we obtain from (|44ll 

and the Matsubara sum (|51|) reduces to an integral, 

m=— oo J —oo 

yielding 



G ™ sf + n( w a, + 4n2)- 

Correspondingly, Eq. 1)46)) becomes 

" m= — oo J oo 

where 



= 0. 



denotes another auxiliary frequency. 
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C. Loop Order I = 2 



The higher loop orders of the effective potential with / > 2 consist of all one-particle irreducible vacuum diagrams 
with the propagator 



G{t 1 ,t 2 ) 



defined by the identity 



and the vertices 

2 



rP i 

J dr 2 G-\n,r 2 ) G(r 2 , ts) = -p *(ti - r 3 ) 



3 = — 



r(i rP rP r/3 

D I dn dr 2 / dr 3 ... dr n 
Jo Jo Jo Jo 



S n A[a] 



5a(Ti)Sa(T2)Sa(T 3 ) . . . 5a(T n ) 
For instance, the two-loop contribution is given by the Feynman diagrams 



(57) 



(58) 



(59) 



'{t)=c?o 



K ( ff 2 Vo) = ioo 



i 



i 



8 — ' 12 - - m 

In order to evaluate (|6U[) . we need the third and fourth functional derivatives of the Euclidean action l|41[l . They are 
given by 

6 3 A\a] 



and 



8a(Ti)Sa(T 2 )Sa(T 3 ) 
S 4 A[a] 



- 64z Gn (n , t 2 ) Gn (t 2 , t 3 ) G n (r 3 , t x ) 



(61) 



<t(t)=(7o 



(5(7 ( Tl ) £ (T ( T 2 ) S a ( T 3 ) J (7 ( T4 ) 



= -256 



c(i")=er 



Gq(ti, t 2 )Gu{t 2i T 3 )Gn(T 3 , T4)Gn(T 4 , ri) 



+G ? n(ri,r 2 )G ! n(T2,r4)G ! n(T4 ) r3)Gn(T3,ri) + Gn(n, t 3 )Go(t 3 , t- 2 )Gq(t 2 , t 4 )Gq(t 4 , ti) 

where the explicit form of the correlation function Go (ti , t 2 ) at zero temperature follows from l|49|) and l|53[l : 

1 



Gn(ri,r 2 ) 



-Q\Ti — T2 



20 



(62) 



(63) 



Furthermore, we have to solve l|58[) for the propagator l|57|l . Performing at arbitrary temperature the decomposition 



1 oo 

G(n,r 2 ) = - VJ G„ 



-iw m (Ti-r 2 ) 



m— — oo 



the coefficient G„, turns out to be 



G,, 



1 



£> Grn' ' 

Using H54JI and (|65|l . we evaluate the Matsubara sum l|64(l at zero temperature according to l|53|l and obtain 



G(n,r 2 ) = — 



'Cn-*)-^e-*~l 



(64) 



(65) 



(66) 



From l|59[) we read off that each vertex is of order D, whereas each propagator is of order 1/D due to (JBfjJ. Thus 
both Feynman diagrams in H6U|) are, indeed, of order 1/D. The first and second Feynman diagram in l|6U|) lead to the 
expressions 



Kg Wo) 



1 



.9 2 



/3£> [ 20 4 
2.9' 



2/ 2 (20) + / 2 (40) 



2.9 



5 



2/3(0, 0,0 + 0)+ I 3 (2Q, 20, 0) 



6 2 
1 

Jd 



2J 4 (0 + 0, 0, 0, 0, 0, + 0) + J 4 (0, 0, 0, 0, 0, 0) 

-1- 7 3 (20, 20, 20) - -^L J 4 (0, 0, 0, 20, 0, 0) + -^L- J 5 (0, 0, 0, 0, 0, 0, O, 0, 0, 0) 
30 6 iK ' ' ' 7 8 2 >'>>>'»'« ; 



16.9 



30 9 3 



=- J 6 (0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) 



(67) 



(68) 
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respectively. Here we have introduced an abbreviation for multiple integrals with respect to imaginary times: 

rP P0 rP inn \ 

I„ (fti 2 , . . . , fti„, ^23, • • • j &2n, ■■•)=/ dr i / dr 2 . . . dr n exp - V] V] Ojj \n -Tj\\ . 

Jo Jo Jo y i=ii= i+ i / 

Considering the low-temperature limit (3 — > oo, the first three of these integrals read in closed form |29|: 



"12 



4/?(fi 12 + n X3 + ft 2 ; 



2:Sj 



j 3 (ft 12 , n 13l ft 23 ) - (fii2 + nia)(nia + fi23)(fii3 + ^ 

^4(^12, ftl3, ftl4, ^23, ^24, ^34) = 2/3 



1 



ftl2 + ftl3 + ^24 + ^34 
1 

S2i2 + "14 + ft 2 3 + ^34 
1 

ftl3 + fti4 + ft 2 3 + ^24 



1 



1 



ftl2 + ftl3 + ftl4 f2l4 + f^24 + ^34 
1 1 



1 



1 



"12 + ^23 + ^24 "14 + ^24 + ^34 
1 1 



ftl2 + ft23 + ^24 ftl3 + ^23 + ^34 
1 1 



ftl3 + ^23 + ^34 ftl4 + ^24 + ^34 



ftl2 + ftl3 + ftl4 ftl3 + ^23 + ^34 
1 1 



ftl2 + ftl3 + f^i4 f2i2 + ^23 + "2 



(69) 



(70) 
(71) 



(72) 



Furthermore, in order to evaluate l|68[l . we need one particular integral with respect to five and six imaginary times, 
respectively: 



/ 5 (ft,ft,ft,o,ft,o,ft,ft,ft,ft) 



/?(3ft 5 + 42ft 4 ft + 227ft 3 ft 2 + 568ft 2 ft 3 + 656S10 4 + 288ft 5 ) 



i 6 (Q, ft, n, o, o, ft, o, ft, o, o, o, ft, n, n, n) 



2ft 2 (ft + ft) 2 (ft + 2ft) 4 (ft + 4ft) 
3f3(n 5 + un 4 n + 7sn 3 n 2 + ieoft 2 ft 3 + i36ftft 4 + 32ft 5 ) 



nn 2 (n + n) 2 (n + 2ft) 4 (ft + An) 

From (|44|l and 167|) - (|74|l we read off the effective potential for zero temperature up to the order 1/D: 



(73) 
(74) 



KffM = D ( -j . 2 



al n\ n 



n + — 

2 D 



5g 2 g 3 (ft 3 + 4ft 2 ft + 44ftft 2 + 128ft 3 ) 
4ft 5 + 



4ft 8 ft(ft + 2ft) 2 
g 4 (3ft 5 + 31ft 4 ft + 150ft 3 ft 2 + 392ft 2 ft 3 + 656ftft 4 + 480ft 5 ) 

n 9 n 2 (n + n)(n + 2ft) 3 (ft + An) 

4g 5 (3ft 5 + 42ft 4 ft + 227ft 3 ft 2 + 568ft 2 ft 3 + 656ftft 4 + 288ft 5 ) 
ft 10 ft 2 (ft + ft) 2 (ft + 2ft) 4 (ft + 4ft) 

16<7 6 (ft 5 + 14ft 4 ft + 73ft 3 ft 2 + 160ft 2 ft 3 + 136ftft 4 + 32ft 5 ) 
ft n ft 4 (ft + ft) 2 (ft + 2ft) 4 (ft + 4ft) 



(75) 



The ground-state energy of the anharmonic oscillator J5J is found by extremizing the effective potential l|75|l with 
respect to the background ctq while taking into account the auxiliary frequencies i|45|) and 156fl . 



D. Weak-Coupling 

In order to cross-check our large- D result, we specialize now to the weak-coupling regime where the extremal 
background is expanded according to 

(T = —i (sig + s 2 g 2 + s 3 g 3 + s A g A + s 5 g 5 + . . .) . (76) 



10 



Inserting l|76|) into the vanishing first derivative of l|75(l and re-expanding in g, we obtain a system of equations which 
are solved by 

si = h> (77) 

S2 = -^-^D' (78) 
5 45 25 

4 59 73 

54 _ - 2^D~ZA°D^ (80) 

231 19503 9823 275 

55 ~ TaJis + 128^13L» + 16cj 13 D 2 + 2w 13 L> 3 ' ' ' ' ' ^ ' 

Inserting (|76l) - (|81[) into (|75l) and re-expanding in g = Dg yields again the weak-coupling series (|11|) . where the 
weak-coupling coefficients tk of Tableware reproduced in the first three terms of order D k+1 , D k , and D k ~ 1 . 

E. Strong-Coupling 

We now derive the large- D expansion for the ground-state energy in the strong-coupling regime. There, we use for 
the extremal background the ansatz 

<ro = -ig 2/3 (Si + S 2 g- 2/3 + ■■) (82) 
and find for the leading coefficient the expansion 

« - 5^ + (*""■ S 1 " - *") ± + (« ■ S" 2 - ?) " 3 + O (±) . (83, 

From (|75|l , l|82|) , and l|83|l we then obtain the strong-coupling series (|17|) where the leading strong-coupling coefficient 
has the large- D expansion 

oo 

b Q = J2B k D^- k , (84) 



fc=0 



with 



3 • 2 1 / 3 

So = — — , (85) 
31/2 _ 0V2 

Bl = T/6 ~ 0.2831607943221791188446646047948820365123, (86) 

239 

S 2 = ——. — ^ « -0.1537760559399284913195761085499705701590. (87) 

18 -22/3 (25 + 12 -e 1 ^) v ; 

Figure |21 shows the leading strong-coupling coefficient l|84|) - (|86|l up to the first two orders plotted as a function of the 
dimension D. 



VI. LARGE-D EXPANSION FROM VPT 



In this section, we study how the large-D expansion (|84|l follows from VPT. Thereby, we numerically determine 
the large- D expansion coefficients up to B§ with high precision. 
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b (D) 




FIG. 2: Strong-coupling coefficient bo(D) versus dimension D. Solid and dashed lines represent leading and subleading results of 
the large-D expansion (I84H - 186I I. respectively. For D = 1, the dot represents the earlier result from Refs. 0,0; for D — 2, 3, 10, 
the dots indicate the highly accurate VPT values obtained from the 80th order, as given in Table [TTT1 



Coefficients Bi and B2 



To first order, VPT gives l|28|l. and we obtain the expansion (|84|l simply by expanding the cubic root in 1)28(1 in 
powers of 1/D, The first three coefficients in the expansion l(84|l read 



R (i) _ 
B - 



3 ■ 2V3 



2 l/3 



2 l/3 



(88) 



The leading coefficient ((8511 is reproduced exactly while the next two subleading coefficients B\ and B2 are missed 
by 11.2 % and 36.6 %, respectively. In the second order of VPT, Eq. ilol'l) yields 



3L>f7 2) 
16 - 



D(2 + D) D{10 + 9D + 2D 2 ) 



2R' 



(2)> 



8(^' 



(2) 



(89) 



which has to be extremized with respect to the variational parameter fig . Setting the derivative of 1(89() to zero, the 
large- D expansion 



n 



(2) 



D l/3 C (0,2) + 



(0,2) 



a 



(0,2) 



leads to a system of equations for the expansion coefficients whose solutions read 



^(2) _ 9 l/3 
°0 — z ' 



(2) _ 13-2V3 
Cl "^2 _ 



> = - 



113 • 2 1 / 3 
288 



(90) 



(91) 



Thus, by inserting the optimized variational parameter 1(90(1 . ((91(1 into 1(89(1 and expanding in 1/D we obtain the first 
three coefficients in 1(841) to second order of VPT: 



B, 



(2) 



3-2V3 



B 



(2)_7-2V3 



32 



(2) 



71 ■ 2 1 / 3 
768 



(92) 



(2) 

The leading large- D coefficient B remains the same, whereas the error of the subleading coefficient B\ ' is reduced 

to 2.67 % and the next subleading coefficient B^ now comes out with an error of 24.2 %. Figure |31 shows that the 

VPT approximants b[ N ^ and converge exponentially fast to their exact values B\ and Bi in ((86(1 and ((87(1 . If 

we had not known the exact analytic result for the subleading coefficients B\ and J5 2 , we could have extracted its 
value from the VPT approximants as follows. Figure |21 shows that the VPT approximants for odd and even orders 



converge towards the exact value independently. Extrapolating the values of B\ 
the limiting values 



(odd) 



and B 



(even) 



B\ 
B\ 



(odd) 



0.28316079432217911884466460479488203808, 
i 0.28316079432217911884466460479488203575 . 



for N — ► 00 leads to 

(93) 
(94) 
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FIG. 3: Deviation of the VPT approximants for B[ N ^ and B^ from their exact values Bi and B2 in Eqs. I86H . 1871 . 
respectively. Odd orders are represented by circles; even orders by triangles. The dashed lines represent fits of the data to 
exponential functions. 



Assuming that the exact value lies within this interval, we obtain from this extrapolation method the result 

B (oxtra P ) = 0.2831607943221791188446646047948820369(24) . (95) 

An analogous procedure can be applied to extract a numerical value for the subsequent coefficient B 2 - Applying the 
extrapolation method for odd and even orders, we obtain the result 

B (oxtra P ) = _ . 153776055939928491319576108549961(60) . (96) 



B. Coefficients up to Bg 

Figure 0| shows that the VPT approximants for B\ and B 2 rapidly converge towards their exact values. For the 
large- D expansion coefficients of higher orders in l/D, however, the VPT approximants first fluctuate and then enter 
the regime of exponentially fast convergence. If we want to obtain the coefficients of the expansion (|84l) to higher 
orders and with good accuracy, we must therefore drive our VPT calculation to high orders. To this end, we specialize 
the expression for the leading strong-coupling coefficient as given in l|32|) in such a way that we can read off the 
corresponding large- D expansion. The weak-coupling coefficients for the ground-state energy ek can be expanded in 
powers of the spatial dimension D, 



fc+i 

3 = 1 



< >J 



(97) 



whereas the variational parameter can be expanded in 1 / D 



A I 



(98) 



Here, M denotes the highest order to which we seek to drive the large- D expansion (|84H . Furthermore, recall the 
multinomial expansion 



x 2 



££-£*h£ 



a-k 



n 



ai— a2— 



fe=l 



where n is an integer and where the Kronecker delta 8ij is written as 5(i,j) 
generalized to hold for real exponents and infinite series as follows: 



(l+x 1 +x 2 + ...) c 



£ 

m=0 



oo oo 



££•••£* 

a\— a2— a m — 

r(a- 



111 ' ' 

1=1 



(99) 



This multinomial expansion can be 



F(a - oi - a-x — . . . — a m + l)ai! a 2 ! 



n 



1=1 



(100) 
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FIG. 4: Logarithm of the relative deviation of the VPT results Bf_ ' for the coefficients Bi through Bg in the large- D expansion 
1841 versus the order N of VPT. The lowest curve (circles) is for B\ N , the uppermost curve (squares) is for Bg , intermediate 
curves are for B^ through Bg N ^ (i?2 ' : triangles, B% : diamonds, B^ N ^: pentagons, B^ N ^ : upside-down triangles). For B^ 
and B^ the exact values from (1861 and 187H were used in order to determine the relative deviation. For B% through Bg N ^ 
we used our extrapolation results from Table iTVl 



Using H97|l and (|98|l in (|32(1 and applying the multinomial expansion (|l()()(l . we obtain 



D 



A I 



j=o 



fe+1 



(0,N) 



N 



D~ 3 



where the coefficients are given by 



fc=i 



N -k 



M+k Mm{j-1, L(M+i-fc)/2J-l} 
j=l j=2 m=l 



,( fe ) K ( N )f)j-2m 
j—rn m.k 



(101) 



K 



(AO 



AI M 
ai— a2— 



M 



(3k 



- ai + a 2 



(3fc-2)!ai!a 2 ! 



n 



(0,7V)' 



(0,JV) 



(102) 



The summation boundaries in p()l|l are determined by the two conditions that i) 



we can neglect contributions containing coefficients C^' N ^ 
derivative of the leading strong-coupling coefficient with respect to the variational parameter d / d 



- for j < 1 V j > k + 1 and ii) 
with k > M . Note that a similar expansion holds for the 

Using 

(jlOip . I|102|) we can efficiently calculate the VPT approximants to high orders. In each order of VPT, the leading 

large-D coefficient Bq comes out exactly. Furthermore, the expansion coefficients C^ '^ for the variational parameter 
are found by solving linear equations. Figure 0] shows the convergence behavior of our VPT approximants for B\ 
through Bq up to N — 100. By extrapolating the VPT results for k > 3 in the limit N — > oo, we are able to 

determine the coefficients in the large- D expansion to high accuracy. The numerical results are shown in Table ITVl 



VII. SUMMARY 



We have determined the weak-coupling series of the ground-state energy for the /^-dimensional anharmonic oscillator 
P)l and used variational perturbation theory to extract the coefficients of the large- 1? expansion to higher orders than 
accessible by standard field-theoretic methods based on the Hubbard-Stratonovich transformation. 
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B 3 


0.1098507701648367991179224418(38) 


B 4 


-0.02886373198599697649759(11) 


B 5 


-0.1695976321261828993(92) 


B 6 


0.51902738902696(86) 



TABLE IV: VPT results for the coefficients in the large- D expansion 1841 of the leading strong-coupling coefficient bo. 
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